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Abstract 

Fractional diffusion equations for three-dimensional lattice models based on fractional- 
order differences of the Griinwald-Letnikov type are suggested. These lattice fractional 
diffusion equations contain difference operators that describe long-range jumps from one 
lattice site to other. In continuum limit, the suggested lattice diffusion equations with 
non-integer order differences give the diffusion equations with the Griinwald-Letnikov 
fractional derivatives for continuum. We propose a consistent derivation of the fractional 
diffusion equation with the fractional derivatives of Griinwald-Letnikov type. The sug¬ 
gested lattice diffusion equations can be considered as a new microstructural basis of 
space-fractional diffusion in nonlocal media. 

PACS: 05.20.-y; 45.10.Hj; 61.50.Ah; 02.70.Bf; 


1 Introduction 

Fractional calculus and differential equation of non-integer orders P 121 El SI E] have a long 
history that is connected with the names of famous scientists such as Liouville, Ricmann, 
Griinwald, Letnikov, Marchaud and Riesz. Derivatives and integrals of non-integer orders have 
a lot of application in different areas of physics pm ei ei unj. Fractional calculus is a powerful 
tool to describe processes in continuously distributed media with nonlocal properties. As it was 
shown in [mug, the continuum equations with fractional derivatives are directly connected 
[7] to lattice models with long-range interactions. The lattice equations for fractional nonlocal 
media and the correspondent continuum equations have been considered recently in [26, [571 
[28] . Fractional-order differences and the correspondent derivatives have been first proposed 
by Griinwald [T9] and by Letnikov [50] ■ At the present time these generalized differences and 
derivatives are called the Griinwald-Letnikov fractional differences and derivatives PEI p El]. 
One-dimensional lattice models with long-range interactions of the Griinwald-Letnikov type and 
the correspondent fractional differential and integral continuum equations have been suggested 
in [23j . The suggested form of long-range interaction is based on the form of the left-sided and 
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right-sided Griinwald-Letnikov fractional differences. A possible form of lattice vectors calculus 
based on the fractional-order differences of the Griinwald-Letnikov type has been suggested 
in [22] , In this paper, we apply this approach to describe diffusion on lattices with long- 
range jumps and to derive fractional diffusion equations for nonlocal continuum with power-law 
nonlocality. 

The diffusion equations describe the change of probability of a random function in space 
and time in transport processes, and it usually have the form of second-order partial differential 
equation of parabolic type. Unfortunately for complex nonlocal media, the usual second-order 
diffusion equation cannot adequate describe real processes. For example, the diffusion processes 
with the Poissonian waiting time and the Levy distribution for the jump length cannot be 
described by equation with second-order derivatives with respect to coordinates. The Levy 
distribution describes the Levy flights [2U 25] that are random walks, where the jump lengths 
have probability distributions with heavy-tails. The Levy motion can be described by equation 
with spatial derivatives of non-integer orders fi, [25]. In this case, the moment of order <5 for the 
Levy motion has the form (|x(t)|' s ) ~ t s ^, where 0 < 5 < /i < 2. Usually the space-fractional 
diffusion equations are obtained from the second-order differential equations by replacing the 
second-order space derivatives by fractional-order derivatives. Fractional diffusion equations 
with coordinate derivatives of non-integer order have been suggested in [II]. The solutions 
and properties of these equations are considered in [15] XT] [13]. The diffusion equations with 
fractional coordinate derivatives was also considered in [8] [9, 18] . 

In this paper we propose a consistent derivation of the space-fractional diffusion equation 
with Griinwald-Letnikov derivatives of non-integer orders by using lattice models with long- 
range jumps that is considered as new microstructural basis to describe fractional diffusion 
processes in nonlocal media. In this paper, we suggest a lattice equation for probability density 
of particle in unbounded homogeneous three-dimensional lattice with long-range jumps between 
lattice sites. We prove that continuous limit for the suggested lattice diffusion equation gives the 
space-fractional diffusion equation for non-local continuum. The fractional differential equation 
for continuum contains the Griinwald-Letnikov derivatives on non-integer orders. Continuum 
models can be considered as a continuous limit of lattice models, where the length-scales of a 
continuum element are much larger than the distances between the lattice sites. The suggested 
approach to derive the space fractional diffusion equations can serve as a microstructural basis 
to describe the spatial-fractional diffusion processes. 


2 Fractional diffusion equation for lattice 

The differences of fractional order and the correspondent fractional derivatives have been in¬ 
troduced by Griinwald in 1867 and independently by Letnikov in 1868. 

Differences of non-integer orders are defined as a generalization of the integer-order difference 

V™ ± /0) = ^—t- 7 ^—t f(xTna), (a G R+, m G N). (1) 

’ Z —' n\ (m — n)\ 

n= 1 x 7 
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The Fourier transformation F of the fractional-order differences of Grunwald-Letnikov type (see 
Section 20 of [H Ej) lias the form 

F {V“, ± /(U} (k) = (1 - e ± “ t )“ T{f] (fc). (2) 

The differences of integer orders ([T]) are defined by the finite series. The differences of non¬ 
integer order a G R + are defined as infinite series (see Section 20 of mm- Fractional-order 
differences of Grunwald-Letnikov type are defined by the equation 


v“. ± /M = E 


-l) n r(a + l) 


r(n + l)T(a -n + 1) 

n=0 \ \ / 


f{x =F na), (a > 0). 


(3) 


The difference V“ + is called left-sided fractional difference, and is called the right-sided 
fractional difference. We note that the series in (j3J) converges absolutely and uniformly for 
every bounded function f(x) and a > 0. 

Using the fractional-order differences ()3l) , we can consider the derivatives of non-integer 
orders. The left- and right-sided Grunwald-Letnikov derivatives of fractional order a > 0 are 
defined by the equation 


GL 


Dl±f(x) 


lim 

CL —^0-b 


v; ± /m 


(4) 


For integer values of a = m G N the Grunwald-Letnikov derivatives (j3J) are equal to the usual 
integer-order derivatives up to the sign in the form 


gl DIJ(x) = (±1)” (5) 

We can note that the Grunwald-Letnikov fractional derivatives coincide with the Marchaud 
fractional derivatives for the functions from L r (R), where 1 ^ r < oo (see Theorem 20.4 in 

mm- 

Let us consider three-dimensional unbounded and bounded lattices. Physical lattices are 
characterized by space periodicity. For unbounded lattices we can use three non-coplanar 
vectors a 1; a 2 , a 2 , that are the shortest vectors by which a lattice can be displaced and be 
brought back into itself. Sites of this lattice can be characterized by the number vector n = 
( 77 . 1 , 77 . 2 , 713 ), where rij ( j = 1,2,3) are integer. For simplification, we consider a lattice with 
mutually perpendicular primitive lattice vectors slj, (j = 1,2,3). We choose directions of the 
axes of the Cartesian coordinate system coincide with the vector a r In this case a j = , 

where a 3 = |a^ | >0 and e 3 are the basis vectors of the Cartesian coordinate system. This means 
that we use a primitive orthorhombic Bravais lattice. Then the vector n can be represented as 
n = riiei + n 2 e 2 + n 3 e 3 . 

Choosing a coordinate origin at one of the lattice sites, then the positions of all other site 
with n = (ni, n 2 , 7i 3 ) is described by the vector r(n) = riiai +n 2 a 2 + ?i 3 a 3 . The lattice sites are 
numbered by n, so that the vector n can be considered as a number vector of the corresponding 
particle. We assume that the positions of particles in the lattice coincide with the lattice sites. 
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The distribution function, which describes probability density for lattice site n, will be denoted 
by /(n, t) = f(ni,n 2 ,ri 3 ,t). This function satisfies the conditions 

+oo +oo +oo 

E E E f{n i,n 2 ,ra 3 ,t) = 1, f(ni, n 2 , ra 3 , t) > 0 (t > 0). (6) 

m=—oo n2=—oo ns=—oo 


To describe dynamics of the distribution function f(n,t) in the lattice models with long- 
range jumps between sites, we define fractional-order difference operators of the Griinwald- 
Letnikov type in the direction e, = /1 a■/ ( of the lattice. Fractional-order difference operators 


of the Griinwald-Letnikov type for unbounded lattice are the operators GL D^ 
the function /(m, t) as 



that act on 


gl O± 





-y +00 

— 'V. GLr % K - m i) /( m > t) ( a j > °> 3 = !, 2 > 3), 

a.; 

J rrij=— oo 


(7) 


where the kernels GL K^{n) are defined by the equation 


gl k ± (n] = (-!)"T(1 + gj) (H[n\ ± H[-n\) 
° j 2T{\n\ + l)T(l + a j - \n\) 


('OLj >0, n 6 Z), 


( 8 ) 


and T(E) is the gamma function, H[n] is the discrete variable Heaviside step function that is 
defined as H[n\ — 1 for n > 0, and H[n] =0 for n < 0, where n6Z. The parameter a 3 is called 
the order of the operator. It should be notes that the definition of Z/’[0] = 1 for discrete variable 
Heaviside function is significant for us, since it allows us to write the kernels GL K^(n) in the 
simple form without allocating repeated zero terms. Fractional-order difference operators (171) 
can be called a lattice fractional partial derivative in the direction ej = /1 a^ |. 

It should be noted that one-dimensional lattice models with the long-range interaction of the 
form GL K+{n) and correspondent fractional nonlocal continuum models have been suggested 
in p3]J. The lattice operators (JTJ) recently have been proposed in [22] . 

It is easy to see that the kernels GL K+{n) and GL K~(n) are even and odd functions such 
that GL K^(—n) = ± GL Jt^(?r). The form of the lattice operators (JTJ) can be defined by the 
addition and subtraction of the Griinwald-Letnikov fractional differences 


Gi D± 


a 


/(m,t) = 


E 


(-l) ro * r(l + a) 


777.0 =0 


= n 2 r (m j + 1 )r( 1 + a-m j ) 


/(n — rrijej,t) ± /(n + mjej , t) 


(9) 


We should note that in equation (J^J) the summation is realized over non-negative values rrij, in 
contrast to the sum over all integer values in equation (JT]) . 

For bounded physical lattice models the fractional-order difference operators also can be de¬ 
fined. Fractional-order difference operators of the Griinwald-Letnikov type for bounded lattice 
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with m' < nrij < m 2 are the operators that act on the function /(m, t) as 


GLpdz 


Gtj 


L 3 


=^ X alK %( n i(j = 1,2,3), 


( 10 ) 


7 1 

rr? ■ — m 


where the kernels GL K±.{n) are dehned by equations (JSJ). The suggested forms of fractional 
difference operators for bounded physical lattice models are based on the Griinwald-Letnikov 
fractional differences on finite intervals (see Section 20.4 of mi). For the finite interval [xj, x'j ], 
the integer values m], m 2 and rrij are dehned by the equations 


m) = 



2 

<M -c- 

1 _ 

, mj = 

1 

X 

1 _ 


, m j = 

— 

— 

, a j . 

. a j. 


- a j. 


( 11 ) 


where the brackets [ ] of (1X11) mean the floor function that maps a real number to the largest 
previous integer number. 

Using the semigroup property for fractional differences of non-negative orders (see Property 
2.29 in 0). it is easy to show that the semi-group property holds for the fractional operators 
(JTj) in the form 


Gi O^ 


Gi O± 


= Gi O± 


OL -f- f3 
J 


(a > 0, (3 > 0). 


( 12 ) 


Using this equation, it is easily to prove the commutativity and the associativity of the lat¬ 
tice operator (JTJ) of the Griinwald-Letnikov type. The commutativity and associativity of the 
fractional operators (j7J) of the Griinwald-Letnikov type for different directions are obvious. 

To describe isotropic physical lattices we should use the difference operators GL 1D)^ and 

^ L 1ED^ with orders aq = a for all j = 1, 2, 3. 

Let us give possible equations for distribution function /(n, t) on unbounded and bounded 
lattices. For unbounded homogeneous lattice the diffusion equation for probability density can 
be considered in the form 


df ( n ; t) 
dt 


X: 

i —1 


gl b ± 


/(m,t),+ 9i i 

i,j =1 


cl d ±,± 


oti (3j 

i j 




(13) 


For bounded lattice we should use the fractional difference operators (TTOj) . and the correspondent 
analog of equation (fT3l) has the form 


df(n,t) 

dt 


Gim\± 


i =1 


Oii 

i 


/( m X)i+ 9i i 

*,j=l 


gl D ±1± 


OiiPj 

i J 


/(m,t). 


(14) 


Equations m and (Hi are the three-dimensional lattice diffusion equations that describe 
fractional diffusion processes with the lattice jumps. Here /(n, t) is the probability density 
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function to find the test particle at site n at time t. The italics i,j € {1; 2; 3} are the coordinate 
indices, g l and are lattice coupling constants. 

The first and second terms of the right hand side of equations (1T31) and (Tl4|) describe the 
particle drift and diffusion on the lattice. These correspondent kernels describe the long-range 
drift and diffusion to n-site from all other m-sites. The parameters cq and /3j in the kernels 
are positive real numbers that characterize how quickly the intensity of the drift and diffusion 
processes in the lattice decrease with increasing the value n — m. The kernels K± (rij — rrij), 
where j = 1, 2, 3, describe long-range jumps in the direction with lattice step length | rij — m 3 \ 
in the lattice. In equation (fT3j) . we use the combination of the lattice operators 


GLp±,± 

ca ■ 

G 

= GL D ± 

Oii 

GLp± 

■ 

-1 


J 


. i . 


L J J 


(15) 


where i, j take values from the set {1; 2; 3}. The action of the operator (fl5l) on the lattice 
probability density /(m, t) is 


GLp±,± 


OLi Pj 

i J 


+00 +CO 


/(m, t)= K ti( n i - m i ) K tMo - m i) /( m > *)■ ( 16 ) 


rm=—c 50 rrij =—oo 


Equations (1T31) and (H4li describe fractional diffusion processes on the physical lattices, where 
long-range jumps can be realized. The lattice diffusion equations (TT3l) and (TT4|) can be consid¬ 
ered as lattice analogs of the fractional diffusion equations for the processes with the Poissonian 
waiting time and the Levy distribution for the jump length. 

Suggested lattice equations (TT51) and (fTTD can be considered as master equations that allow us 
to describe time-evolution of particles and quasi-particles on lattice since evolution is modelled 
as being in exactly one of countable number of lattice sites at any given time, and where 
switching between sites is treated probabilistically. These equations are differential equations 
for the variation over time of the probabilities that the particle occupies each of the lattice 
sites. 


3 Fractional diffusion equation for continuum 


To describe fractional diffusion in the nonlocal continua, we should use fractional derivatives 
with respect to space coordinates instead of the lattice operators. Continuum analogs of the 
fractional-order difference operators of the Griinwald-Letnikov type are the fractional deriva¬ 
tives of Griinwald-Letnikov type. 


Fractional-order difference operators Gi O^ defined by (0 are transformed by the con¬ 
tinuous limit operation into the fractional derivative of Griinwald-Letnikov type with respect 


to coordinate Xj in the form 


lim 

CLj —^OH - 


GL H)f 


CH-j 


L J 



gl B± 


a.j 


L J 




(17) 
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where Gi O^ 
are defined by 


are the continuum fractional derivatives of the Griinwald-Letnikov type that 

GLj^a 


gl b± 


GL T^a 


± G L D 


(18) 


which contain the Griinwald-Letnikov fractional derivatives GL D “ ± with respect to space co- 

(- 1 )" 1 ' r(a + 1 ) 


ordinate x 3 that can be written as 


GL na 


D“ >± /(r,t)= lim 


E 


a ^°+ M“ “ 0 r ( m i + 1 ) r ( a “ m i + X ) 


/(rT^-ajd), (a > 0). (19) 


This statement can be proved by analogy with the proof for lattice model with long-range 
interaction of the Griinwals-Letnikov type suggested in [23] , 

It is important to note that the Griinwald-Letnikov fractional derivatives coincide with the 
Marchaud fractional derivatives (see Section 20.3 in BED for the functions from the space 
L r (R), where 1 ^ r < oo (see Theorem 20.4 in M). Moreover both the Griinwald-Letnikov 
and Marchaud derivatives have the same domain of definition. The Marchaud fractional deriva¬ 
tive is defined by the equation 

i r oo Ai>±f( r t) 

—^ TI - dzj, (0 < a < s), (20) 


M D a x ?f{v,t) = 


a(a, s ) 


where i s the finite difference of integer order s, 




-1 ) k s\ 


k =0 


k)\k\ 


/(r — k Zjej, t), 


and a(a,s ) is 


a(a, s ) = 


^ [G i-0 


s-i 


-A- 


( 21 ) 


( 22 ) 


a Jo (ln(l/£))< 

Using (JSJ), we can note that the derivatives (fTHj) for integer orders a = n E N have the forms 

iil-D-a" (23) 


g1 d± 


dxj 


Therefore the continuum fractional derivatives GL f 


'c 


are the usual derivatives of integer 


order n for even values a only, and the continuum operators Gi D^ 
integer order n for odd values a only. 

For bounden lattices, the fractional-order difference operators 
are transformed by the continuous limit 


lim 

CLj -^0“|- 


gl D ± 

B % 


a 
L 3 


3 /(m,t) = Gi O± ' 


a 3 

L 3 


are the derivatives of 
defined by (TTOl) 

(24) 
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into the continuum fractional derivatives of the Griinwald-Letnikov type with respect to space 
coordinate x 3 . 


b L ®c 


GL n a 
x J 


± % L D a 


Xj Xj,— I ’ 


(25) 


which contain the Griinwald-Letnikov fractional operators defined on the finite interval [xj , xj] 
with xj = mjdj and xj = rrijaj, in the form 


< B LD 2 j ,±fi r ,t)= lim 


M f 

E 


-l) m * r(a + 1 ) 


a j^o+ \ aj \ a ^ 3 r(m j + l)r(a - rrij + 1) 


/(r T mj&j,t), 


(26) 


where 


M+ = 


Xj Xj 


M~ = 


Xj Xj 


(27) 


The suggested forms of continuum fractional derivatives of the Griinwald-Letnikov type allow 
us to consider diffusion processes on bounded areas of nonlocal continuum. 

The lattice diffusion equation (1T3|) in the continuum limit gives the fractional diffusion 
equation with derivatives of non-integer orders with respect to space coordinates. This space- 
fractional diffusion equation for the probability density /(r, t) has the form 


g/M) 

dt 


-£ D ‘(“) 


G1 B5 


2=1 


OLi 

L i 


3 3 




GL n ~- 


i= 1 j =1 


Oii /3j 
i J 


/M), (28) 


where Dija) is the drift vector and Dij(a,{3) is the diffusion tensor for the continuum that are 
defined by the lattice coupling constants g t and g tJ by the relations D^a) = g j, Dij(a, (3) = 2 g tJ . 
Similarly, the diffusion equation for bounded lattice gives the fractional diffusion equation for 
bounded region of continua 


dt 


E>,i 


a 


Gi m)- 

B 


i= 1 


Oii 

L i 


3 3 


f(r,t) + -22E D ‘i( a ’W “ D c~ 


i= 1 j =1 


Oii f3j 
i J 


/( r > t). (29) 


It should be noted that coincidence of orders of fractional derivatives in the first and second 
terms allows us to represent the fractional diffusion equation (I25|) in the form of the space- 
fractional continuity equation 


df{r,t) 


dt 


£ 

2=1 


gl Bc 


Oii 




(30) 


where J* is the probability flow 


\ E D o-( a > £) GL]D> c ^ 

j =i L 3 


/( r > t). 


Ji(r,t) = Di(a ) /(r, t) 


(31) 



























Equation (13T1) can be considered as the fractional phenomenological Fick’s first law for nonlocal 
media. If cq = 1 for all i = 1,2,3, the continuity equation (1301) can be represented as the 
standard form of the well-known continuity equation 

<9/(r, t) = A dJj(r,t) 

dt ^ dxi 

1=1 

where Jj(r,t) is defined by (T5TT) with f3j ^ 1 in general. 

For one-dimensional case with D^a) = 0 and /(r, t) = f(x,t), 
sented in the form 

where K{/j) is the generalized diffusion constant, 

K(») = ±Dn(a,P), (34) 


equation (1251) can be repre- 


(33) 


and V M is the fractional derivative of order /i, 





'CKi" 

- X - 





Oi\ -\- fil 
X 


[i — oi\ -\- f3\. 


(35) 


Here we use the semi-group property of fractional derivatives of the Griinwald-Letnikov type. 
Equation (133|) describes the fractional diffusion processes with the Poissonian waiting time and 
the Levy distribution for the jump length (see Section 3.5 of |25j). In [25j the space-fractional 
diffusion equation (1331) contains the Weyl fractional derivative of order /i, of one-dimensional 
case. The solution of equation (l33j) can be obtained analytically by using the Fox function PL) 2 
(for details see Section 3.5 of [25j and |29j). The exact calculation of fractional moments 25) 
gives 


<M*)I S > 


2(A'M) i -'>T(-y / u)r(i + ^) 
/xr(—5/2) r(i + 5/2) 


(36) 


where 0 < 5 < p. < 2. 

Using equation 02]) , it is possible to demonstrate that the space-fractional diffusion equations 
are connected with continuous time random walk processes with diverging second moment of 
the jump length distribution [[25]. 

If aj = f3j = 1 for all j = 1, 2, 3, then equations (125]) and ([29]) give the well-known second- 
order diffusion equation 


<9/ ( r ) t) = _ n. ( 9/( r ^) , n. *9 2 /(M) 

dt 1 Q x 2 2.^ P Q x Q x 

i= 1 i =1 j=1 J 


(37) 


where D, L = D, ; (l) is the drift vector and D, 3 = ZUj( 1,1) is the diffusion tensor for local 
continuum. 


9 















4 Conclusion 


Lattice analog of the fractional-order differential equations for bounded and unbounded three- 
dimensional lattices with long-range jumps of particles are suggested. These lattice equations 
can be considered as a new microscopic basis to describe the fractional diffusion in nonlocal 
continua. In the lattice diffusion equations, we use the fractional-order difference analogs 
of fractional derivatives, which are represented by kernels that describe long-range jumps of 
lattice particles. The proposed kernels of long-range jumps on the lattice can be considered for 
integer and fractional orders of suggested difference operators. The continuous limits for these 
diffusion equations with fractional-order differences give the continuum fractional derivatives 
of the Griinwald-Letnikov type with respect to space coordinates. The obtained fractional 
diffusion on nonlocal continua can be considered as a continuous limit of the suggested lattice 
diffusion, where the sizes of continuum elements are much larger than the distances between 
sites of the lattice. The main advantage of the suggested approach is a possibility to consider 
fractional-order difference diffusion equations as tools for formulation of a microstructural basic 
model of fractional diffusion in nonlocal continua. The proposed three-dimensional lattice 
diffusion equations can play an important role to formulate discrete models of nonlocal processes 
in microscale and nanoscale. 

It is interesting to generalize the suggested lattice approach to consider lattice Levy flights 
subject to external force fields and the Galilean invariance. This transport process on lattice 
can be described by lattice fractional diffusion equations. We assume that the lattice Levy 
flights in a constant force field are similar to lattice fractional diffusion in a constant velocity 
field by analogy with diffusion processes in continuum models [30] . 

We assume that the proposed lattice approach to the lattice fractional diffusion can gener¬ 
alized to different types of Bravais lattices such as monoclinic, triclinic, hexagonal and rhom- 
bohedral lattices. We also assume that the suggested approach to the fractional diffusion can 
be generalized for lattice models with the fractal spatial dispersion, which are suggested in [31] 
(see also [32, [33]), and the continuum limits of these fractal lattice models can give continuum 
models of fractal media [36] [37] that are described by non-integer-dimensional space approach 
[Ml [35]. 
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